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ABSTRACT 

We propose a simple and solvable mean-field model of the scattering of transverse optic modes by 
Polarized Nano Regions (PNR) in the paraelectric phase of relaxor ferroelectrics. The PNR is assumed to 
be a ferroelectric sphere embedded in the host isotropic medium. The Lagrangian parameters are taken 
to be the same inside and outside the PNR, with exception of the soft-mode gap temperature 
dependence. The interaction of the Transverse (TO) with the longitudinal (LO) optic modes is taken into 
account but the latter is found to be important only at the surface of the PNR. Elementary excitations of 
the system are found to be of two types - Vortex (V) and Quasi Polar (QP). V excitations correspond to 
closed polarization lines or closed TO displacements while QP excitations contain open TO polarization 
lines, with an electric dipole if the net angular momentum 7=1. LO waves are virtually excited only in the 
thin layer near the PNR surface and can be excluded from the boundary conditions without any 
consequence. The final boundary conditions include only V and QP TO wave amplitudes. Dynamical 
equations are solved and a phase diagram is predicted. The phase transition temperature is found to be 
lower for Vortex than for QP excitations. Therefore, the QP condensation (local phase transition) occurs 
before the Vortex one upon cooling. TO scattering by the PNR reveals strong long-wave resonances that 
can be shown to result from shallow localized and quasi-localized states and may be essential to 
understanding the waterfall observed in TO studies by inelastic neutron scattering. 

PACS: 77.80.-e; 77.84.-s; 78.70.Nx; 77.22.Gm 



1. INTRODUCTION 

Over the past several decades, serious progress has been made in understanding the physical 
characteristics of the disordered ferroelectrics so-called relaxors (see for example reviews [1], [2]). 
These technically important materials (PMN, PZN, KTN and other) are characterized by chemical and 
structural disordering. Early on. Burns and Dacol [3] introduced the concept of Polarized Nano Regions 
(PNR) - small ordered (polarized) regions embedded into the host crystal - as the main feature of 
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relaxors. Since then, the validity of this concept has been confirmed by numerous optical, diffuse and 
inelastic neutron scattering studies (see for example [4, 5, 6]). 

An important observation was made by Gehring et al [5, 7] in their study of the Transverse Optic mode 
(soft-mode) dispersion in the relaxor PZN-8%PT in a temperature range where PNR's should certainly 
exist. These authors found that, while perfectly well defined TO peaks were observed in the neutron 
scattering spectra (constant Q-scan) for momenta q>qo, these peaks disappeared for q<qo and the value 
of 1/qo was comparable to the expected PNR size of a few tens of angstroms. In addition, the dispersion 
curves had the shape of a "waterfall". Different explanations of this phenomenon have been proposed 
[7 -r 9]. Gehring et al [8] attributed their results to a sharp step-like increase in the TO damping for q<qo 
(that could be described by an even Fermi function). However, no firm consensus has been reached as 
to the definite origin of the waterfall phenomenon. Another motivation for the present theoretical work 
are our own neutron scattering measurements of the transverse acoustic mode (TA) dispersion and 
damping in the relaxor K^xTa^NbOa (KTN) [10,11]. The shape of the TA peaks exhibit clear deviations 
from a standard Lorentzian form (Q-scans) at several TA momenta. These deviations appear to be more 
or less regular, not easily explained by experimental uncertainty and could be interpreted as evidence 
for resonant scattering of the TA by PNRs. 

It is now generally accepted that polarization correlations (PNR) appear at the so-called Burns 
temperature, Tb, upon cooling. It has also been argued [12] that the correlations appearing at T=Tb are 
dynamic and become static at a significantly lower temperature, T*, corresponding to a local distortion 
or local phase transition. T* is in fact the temperature at which elastic diffuse scattering appears in 
neutron spectra as well as the characteristic relaxor behavior or frequency dispersion of the dielectric 
constant. Consequently, the term polar nanodomains (PND) has been proposed to describe these 
polarization regions. So far, we do not have detail information about the shape of these regions, the 
nature of their boundaries (sharp or smooth), the value of the surface energy etc. In the present work 
(see also [13]), we consider PNRs in the temperature interval T*<T<Tb, but nevertheless suppose that the 
relaxor can be approximated as a two-component system. We then describe theoretically the scattering 
of the transverse optic phonon (TO) by PNRs in the framework of the simplest possible model - a 
spherically symmetric PNR embedded in an isotropic medium. For simplicity, we assume that all 
Hamiltonian (Lagrangian) parameters are the same outside (index 1) and inside (index 2) the PNR, with 
exception of the temperature of the phase transition (Tel, Tc2) and the energy gaps a)o(l), cjOo(2) of the 
TOW soft-mode at zero momentum, q=0. We also assume that all parameters are constant. Mixing of 
transverse and longitudinal waves at the PNR surface is taken into account but we assume the same gap 
parameter for the LO and TO in order to avoid using higher derivatives in our equations. In the 
Lagrangian, we therefore include the interaction of the polarization with the electromagnetic field and 
take into account longitudinal Coulomb forces. However, we describe the system with a critical TO and 
non critical LO behavior because the LO does not usually soften in ferroelectrics and because Coulomb 
forces tend to raise its frequency at small momentum in ionic crystals, leading to a large energy gap Qp 
» cjOo(1), cjOo(2) and only a moderate (non-critical) dependence on temperature. From electrodynamics, 
we find that elementary transverse excitations of the system can be divided into two categories - Vortex 
(V) and quasi polar (QP^). Vortex excitations correspond to closed or loop polarization lines of TO 
displacements and QP excitations contain open TO polarization lines (if the total angular momentum 
j=l). LO waves are virtually excited (Qp» cjOo(1), a}o(2) ) but only in a thin layer near the PNR surface. In 
the limit Qp->°°, we can analytically exclude them from the boundary conditions without restriction on 
the validity of the results. Hence, the final boundary conditions contain only the TO Vortex and QP 
amplitudes. The vortex boundary conditions are standard and correspond to the condition of continuous 
TO displacement and its derivative with respect to the radial coordinate at the PNR surface. Quasi-Polar 



Similar photon excitations are called as magnetic and electric in electrodynamics [16]. 



excitations correspond to the condition of a continuous TO displacement but discontinuous derivative 
with respect to the radial coordinate at the PNR surface. This discontinuity of the derivative results in a 
near-field effect that vanishes for a large size PNR. 

In the following parts of paper we formulate the model and provide its solution for the 
boundary conditions described above. We show that a local instability ("phase transition in PNR") occurs 
at a temperature Ti<Tc2 where the mode frequency goes to zero locally, cjOl=0 (if Tc2>Tci, Tc2 is then the 
critical temperature for the case of a very large size PNR). f^ is found to be lower for the case of a Vortex 
excitation than for a QP excitation. Therefore, the QP condensation (local phase transition) is expected 
to occur before the Vortex condensation upon cooling. Also, f^ is found to be lower for a smaller size 
PNR. For Tc2<Tci, we also find that local TO modes can exist, for which cjOl < cjOo(1), due to reflection by 
the host medium. In the last part of our paper, we calculate the Vortex and QP scattering by the PNR 
and find strong resonance scattering at small TO momentum. All calculations are exact as we do not 
assume that the TO-TO interaction between PNR and host medium is small. A short discussion of the 
results and future perspective are presented in the conclusion section. Some mathematical derivations 
are detailed in the Appendix. 

2. MODEL 

2.1 MODEL DESCRIPTION. 

We consider the dynamics of a spherically symmetric PNR embedded in an isotropic medium. The total 
Lagrangian density L contains the contributions of the TO and LO waves, Lo, the electrostatic field [15], 
Lem, and the interaction between electric field and polarization, /.,„(. 

The total Lagrangian A is written as: 
A^lLd\, L^L^+L^^+L,^, (1) 

L^ ^K^-n^ where /ir„ = 1/2^^(9^/50^ and /7„ = +{AJ 2)dH^)^ +bJ-^^+\I2p^cdI^'^ 

^,.^=1/2(9^,. /axy^+e^ /ex,.), 4=^/yt-i/3^ik^qq (2) 

^..-^(V^)' (3) 

Ant = -^^a'^aV (4) 

Here ^ is the displacement in the optic mode, Po the corresponding density, e' the effective charge. Ao, 
and Bo are parameters describing the optic mode contribution to the potential energy written in a 
similar form to that in elasticity theory [14] and cp is the electric scalar potential. A similar Lagrangian 
was proposed by Hopfield [15] for the description of polariton dynamics. Our Lagrangian is different 
from that proposed by Hopfield in three essential aspects. First, we do not take into account polariton 
effects. Second, we take into account the space dispersion of optical modes (see flo (1) first and second 
terms), and, third, we also take into account the softening of the optic mode in the last term of flo. 

(Oq (I) = al (r-7cj, (Oq (2) = «2 \T-Tc2], Tc^ < Tc2, al < a2 (5) 



A phase transition is expected to occur inside the PNR at a higher temperature than in the host medium, 
Tci< Tc2, resulting in what we called earlier a polar nanodomain. We should also note that Hopfield used 

the expression Lh for the description of the interaction between electric field and polarization, 

^H='^ <^^cK^« =-^^Tnt +«? '^ctr(^cir^) Both expressions Lh and (4) lead to the Same results for 
the case of infinitely large medium. We prefer apply expression (4) which is similar to the standard 
description of interaction between polarization and electric field and polarization [1] and more suitable 
for the boundary conditions formulation. Given the above Lagrangian, the dynamical equations are 
written in the standard form: 

d I dt{5K I S(d^^ (r) / dt)) = SA 1 5^^ (r), ^A / 5(p{r) = (6) 

From Eq.2 we then get: 

PO^'^^a I ^'^ = ^^ * Srad^dH^) + 2Bo * V^ (^^ ) - 

-^O^X -^*^«^ -^oidiv(^)dS^ -2Bo*l i^dS^ (7a) 



and from Eqs.3 and 4: 

-Va(V«^-4^*^a) + l(V^^-4^^*^J^^^ =0 (7b) 

The electrodynamics boundary condition for the tangential component of the electric field is given by: 



V(Pi-n(nV(p^) = V(P2-n(nV(p2) ^^c) 



where dS is an element of the PNR surface and n is a unit vector normal to the surface. The terms 
containing dS (7a, b) lead to particular solutions satisfying the boundary conditions at the PNR interface. 
(7b, c) are the standard boundary conditions of electrodynamics. It is important to emphasize that all 
parameters are the same inside and outside the PNR with the exception of the gaps, coo(l) and coo(2), 
determined by the transition temperatures Td and Td respectively. We now derive separate 
independent solutions inside and outside the PNR in the form of plane waves with frequency oj and 
momentum Ki, /fy- for the optic longitudinal and transverse modes respectively. From Eqs.7a and 7b: 

Longitudinal Optic wave (LOW) 

2 2 2 2 * 

[PgCO -AoKj^ - 4 / 3BoKj^ - p^co^ ]<^£ - iKj^e (p = 0, 

iK^(p -4nQ #L = (^^) 



Solving this system of two coupled equations, we get: 

2 ^ 2 2 2 

CO - AoK , In - 4 / l,BoK , I p - co^ - O. =0, 
L o L o p 

with fip=4;r(e*)2//5o (8b) 

Because the ion plasma frequency, Qp~lS - 75 meV »cjOo, is large, the softening of the LOW can be 
neglected. 

Transverse Optic wave (TOW) 

Yo?-BoKt^ I Po-o^^%T=^ (9) 

For convenience, we write the parameter Aq and Bq in the LOW and TOW dispersion, Eqs. 8b and 9, in 

terms of the speeds C^ Cq. 

Ao^Po[cl-4/3C^ ], Bo^pqC^ , ^^^ ipo (10) 

We also take the density parameter, po, to be equal to 1. 

2.2 MODEL SOLUTION (GENERAL). 

Because the effective potential energies inside or outside the PNR are independent of the coordinates 
(Eqs. 1-4), the solution can be decomposed into independent radial and angular parts. The radial part can 
be written in terms of Bessel (Hankel) functions inside (outside) the PNR with argument rKj and r/C^ for 
the TOW and LOW respectively. The angular part is written in terms of spherical functions. We are 
especially interested in the dynamics of transverse excitations. A similar problem was considered in 
quantum electrodynamics [15]. Our case is similar but does not quite coincide with [15] because the 
longitudinal optic modes are important in our case while they are absent in [15], and because the optic 
mode energy spectrum contains a gap. For convenience, we utilize some results from Refs. [15], [17] and 
[18] concerning spherical vector functions. Because we are interested in scattering by a spherical PNR 
and the total angular momentum is conserved, it is convenient to use an angular momentum technique 
in the following. We describe the optic phonon polarization in terms of the vector spherical functions, 
X,(a), ^=0, ±1, a=l,2,3 (x, y, z): 





1 f'^ 


f^^ 


Z-i 


-72[o'}'0- 


= 



(11) 
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The phonon wave field Yjim contains angular and polarization ("spin") components in the total angular 
momentum y representation: 

J';/^(v,a) = I ^jmlu^lm^^^^u^"^'' ^^^'^ m + ju = M and \ = k/k, 
m,ju ^ 

V =1, / = 7,7 ± 1 if 7 ^ 0; / = 1 if 7 = and M = -j...j (12) 

where k is the momentum, Yjm is a standard spherical function [15-18] with angular momentum / and its 
projection m and d'^imiii is the Clebsch-Gordon coefficient [15]. We should emphasize that we are 
working here in the momentum (not-space) representation. Expression (12) can be written in the form 



of orthogonal normalized spherical vector functions (integration is done over the surface of a sphere of 
unit radius). 

We distinguish several solutions: one longitudinal, with spherical vector Y^'^'jivilv) (superscript -1), 
directed along the momentum k, and two transverse, y'°'jm(v) and y'^'jmM (superscript and +1), with 
different parity). Several useful expressions concerning spherical wave functions [16 -19] are used in the 
paper: 

V ^^^^ ^ V ^^^' V ^"^ ^ ^JJM (^)' V ^"^ ^ '^^JJM ^^>' ^]' 

Yy(vrv=o, Yy(v)=-.— ==^, 

Y.,^ (v) = ^i^^Y. . , ,^ (v) - \- Y. . , ,^ (v), 



(1) 
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7+1 



Y.,^ (v) = ^i^^Y. . , ,^ (v) + ^ Y. . , ,^ (v) (15) 

jM )i27+l 7J+l,Af ' )i27+l 7J-1,A/ ' ' ' 

We note that the functions YjJ'''(v), YjJ^'(v) correspond to Vortex (V) and Quasi polar (QP) excitations 
respectively. A Vortex excitation is a manifold of closed lines of optic mode displacements (similar to the 
magnetic field lines in electrodynamics), and therefore does not possess an electric dipole moment. A 
QP TOW displacement (15) contains a spatial component corresponding to the spherical symmetric 5- 
orbital state but also spin angular momentum, giving a non-zero electric dipole moment after 
integration over angles when j=l. We note also that only the longitudinal wave exists when the total 
angular momentum y'=0. 

Next, in order to fulfill the boundary conditions at the PNR surface, the phonon wave field must 
be expressed in real space. This can be done by means of an expansion of the plane wave in terms of 
spherical harmonics [15]: 

* 2 

exp(/kr) = S giikr)Y (k/k)Y (n), n = r/r, n = 1, 
l,m '^" '^" 

gl{kr)={27t9'^ihi+ii2{kr)l^r. 



J 



Yy7M(k / k )^xpiikr)dnk = g ^ (kr)Y^.^^ (r/r) (16) 



Ji+i/2 is an ordinary Bessel function. As usual, Bessel functions are used in expression (16) inside PNR to 
ensure a finite TOW displacement at the PNR center, r=0; outside the PNR, we should use the first 
(second) order Hankel Hi+i/J^'''^' function for the case of an "outward wave" ("inward wave") instead of a 
Bessel function. According to our earlier notation, the wave vector outside (inside) the PNR is labeled as 
I<=Kt(1) [k=K-!{2}). In the real space representation, the wave fields corresponding to the V°\mM, y'^'jm(v), 
and Y^'^'jm(v) components are transformed into the expressions /^"'jm, T^^^m, ond F^'^'jm, respectively: 



(0) (0) (0) (0) 

^jM = Sjikr)Yjj^(n)^ fj (kr)Yjj^(n), (17a) 



/j^'^^kr)YJj^(n)+/j^'"^\kr)YJ^\n) (17b) 

(-1) j 7 + 1 (-1) 

^^|ff5 [gy-l(^r)-g,-+i(^r)]YJ^(n). 

(-1,-1) (-1) (-1,1) (1) 

fj (kr)Y;^ (")+/) (kr)Yy(n) (17c) 

where we have introduced the coefficients /°''j> f^'^\, /^''^> /^'^> /"''"'''for the expansion of the wave field 
components in terms of the orthogonal basis functions Y'°'jM(n), Y'^'jM(n), Y^'^'jiviin)- 
It is important to note that, when going from the momentum to the real-space representation, a 
transverse vortex (0) excitations converts to a transverse vortex excitation while a transverse quasi-polar 
(+1) converts to a combination of transverse (+1) and longitudinal (-1 ) excitations (17b). Hence, the 
index ±1 in the space representation and ±1 in the momentum representation do not have the same 
physical meaning. Polarization is normally defined with respect to k, i.e. in the momentum 
representation. However, the boundary conditions are applied in real space, so that the conversion from 
momentum to real space representation is necessary. The main physical characteristic of the TOW 
(div(0=0) is of course invariant with respect to the type of representation used. The functions f'°'jm, 
T^*^\m are transverse in the momentum representation (f'°'jm v=0, f'^^^'jm v=0), i.e. with respect to k. In the 
space presentation however, the function f^°\M is still transverse but this time with respect to r and the 
function T^''^\m contains a longitudinal Y^'^^mIT^O term (17b). This term describes near field effects and is 
proportional to 1/r^ in the wave zone (kr»l). At a large distance from the center of the PNR, 

3/2 '— / sm(A:r ) 

kr ^<x, g (kr) « {2n) |^(0 ^^^ (18) 

We refer to /^% as the "Vortex" and to /^% as the "quasi polar", QP, components of the TOW. In a QP 
T^^\m type excitation, both transverse and longitudinal optical mode components will be present as a 
result of applying the boundary conditions (7a), because these BCs are stated for the tangential and 
normal components of the displacements and stresses, with respect to r and not k. By contrast, the 
Vortex f^°\M excitation will correspond to a strictly transverse excitation and will not be mixed with a 
LOW, since it only contains a y°jM term, as seen in expression (17a). The general displacement ft in the 
longitudinal optical wave (8b) can also be written as the gradient of some scalar "potential" Oq (19). 

^, = VCD , with CD ^ G .{r)Y.^ <&,(/)) (19) 
Writing ^ in terms of the potential <^o in Eqs. 7a and 7b, we get: 



2 2 2 ^ 

AO + (o) - 0)^ - n ) / C, O =0 (20a) 
o p L o 

and 

2 2 2 

^f ( r^ ^^ )-^^^Gj ( r )+ ' %-^Gj ( r )=0 (20b) 

y or or y (jj^ 

The LOW displacement (21) contains not only a component xY^'^'j^in) that is perpendicular to surface of 

the PNR but also a component a:y*^'y;^^ (n) that is parallel to it (n=r/r, n^=l). 

dGj dGj (-1) ^,'7(7^ (1) 

E = ^5^nF.,^ + G .V7.,^ = _^ Y.,^ (n) + G . Y.,^ (n) (21) 

^L or jM J jM or jM ^ ' j r }M ^ ' ^ ' 

The near field effects described by the term y*^'jwj ("iG/r in Eq.21 rapidly decay at a large distance r 
from the PNR center ~l/r^. Also, we are interested by the case when the ion plasma frequency in (8b), 
flp, is very large, Dp»(jo. Such a LOW will only be excited within a narrow layer of thickness 6R at the 

PNR surface, SR ~ \/\27rKT (8b). For the case of a relaxor crystal, a typical PNR radius is R~4 -rlS 
l.u. and 6R~1 l.u. if Qp=30meV and Cl=100 meV*l.u. 

We now use the general solution of the model to describe the scattering of an incident 
transverse optic (TO) plane wave by a spherical PNR. Expanding the plane wave in terms of spherical 
harmonics [13], we obtain an expression for the partial wave ct% a'^'j and total, a'"', ct^^^ scattering cross- 
section: 

a ' = ^ (2j+l)sin(^^^'^^)^2;r/yt^ (22) 

j>l J 

where S^- ' and S^- ' are the phase shift for scattering of vortex and quasi-polar excitations 
respectively. 

2.3 BOUNDARY CONDITIONS 

We begin with the simple case of Vortex excitations. The TOW displacement can be written as (17a), 
where/,- is a Bessel function inside the PNR and a Hankel function outside. The last term in Eq. (7a) 
requires that both the optical displacement and its derivative (strain) in the Vortex excitation be 
continuous through the surface 



(0) 
^ = f.{kr)Y.^ (n), ^(r = R - ^r) = ^(r = R + s), e ^ +0 (23a) 

1 d^a ^^P 1 ^fi fj (0) 

We see that the vortex components are not mixed with quasi-polar excitations. 

The vortex excitation boundary conditions can finally be exactly written in a compact form: 

^(5) = 0, ^ = (24a) 

where the auxiliary function MJ(g) represents the difference between the TO amplitudes outside and 
inside the PNR: 

¥(5) = A^^gj^,^ + ^i%- out - A^°\2)g, (24b) 

a'°'(2), A'°'i„, A^°'out3re the Vortex amplitudes respectively inside the PNR, incident (in) and reflected (out) 
from the PNR surface on the outside. 

The derivation of the boundary conditions is more complicated for the case of QP excitations 
(see Appendix). However, these conditions can be written in a compact form when the frequency of the 
longitudinal optical mode is very large. They can be expressed by means of the auxiliary functions (D(g'^') 
and 0(g'""): 

(1) (1) (1) (1) 

<5(G ) = A G^.. +A G„, -A G^,, 

Tlin Tim Tlout Tlout T2 T2 

(-1) (-1) (-1) (-1) 

OfG ) = A^, G^, + A^, G^, - A^^G^^ (25a) 
^ ' Tlin Tlin Tlout Tlout Tl Tl ^ ' 

(1) 

occ? ) = o, 

^ = r ^iG ) (25b) 

CD(G'^')=0 expresses the continuity of the QP TOW displacements at the PNR surface and CD(G*'^'), the 
discontinuity of 0(g'^') derivatives with respect to r at the same surface. This discontinuity is due to the 
excitation of a short-wavelength LO wave, which is important as it leads to the existence of oscillating 
electric charges in a very thin layer near the PNR surface. Expressions for the functions G'^TUn, G'^'jiout, 
G'^'t2; G''^ Tlin, G^'^'jiout^ G'"^'t2 can be found in the Appendix. It should also be noted that G^'^ } describe 
near-field effects and vanishes for large values of the PNR radius. The optic mode displacement, f, is 
described by the relation (Al) in the Appendix. The density of electric charge p=e*div(0 can be 
described by a superposition of spherical harmonics Yjm(^,<p), J^l, so that the total charge inside a thin 

spherical layer is equal to zero, ^pdQ. = 0. Equations (24a, b, 25a, b) allow calculating the phase of 

scattering, its cross-section and all necessary values (for example local mode frequency, space 
distribution of displacement in the Vortex or QP TOW excitations). 



3. RESULTS of CALCULATIONS 

3.1. SOFT-MODE - TEMPERATURE DIAGRAM. 

We first plot in figure 1, the soft mode frequency, a}(T), inside and outside the PNR for the particular 
values of the transition temperatures, T^ and Td, used in the calculations below. General conclusions 
about the TOW dynamics are then obtained in the following subsections. 



200 n 




T(K) 



Fig.l Soft mode energy gap versus temperature for the case of a large radius PNR. Tci=2K, Tc2=125 K, 
al=0.6; o2=l; Co=70. ojoilf=al^{T-Tci}, cjo(2f=a2^(T-Tc2). Numerical values of parameters were 
estimated for the particular physical systems considered and are being used in following calculations. 
Ttransp- Is the temperature at which the soft-mode gaps inside and outside the PNR are the same and the 
soft-mode not scattered by the PNR. Ttransp=(afTcra2^Tc2)/(af-a2^), Ttransp~194.18 K.l- triangle region 
where localized modes are possible; 2 - region where an incident TOW is scattered by PNR; 3 - 
condensed PNR phase. 



We limited ourselves to the temperature region T<T 



transp- 



3.2. SYSTEM INSTABILITY DUE TO VORTEX AND OUASIPOLAR EXCITATIONS 

A PNR embedded in a very large host medium will spontaneously transform to another phase (polar 
nanodomain, PND) through a local phase transition. A large PNR will undergo a phase transition at T=Tc2 
(Fig.l) but a finite size PNR will transform at a lower temperature Ti^st <Tc2. The instability occurs when 
the frequency co of the TOW within the PNR goes to zero. In the following, we consider a temperature 
T>T|nsfyet transformed to a PND. Therefore, we should find the non-zero amplitude solutions (24a, b) 
and (25a, b) corresponding to both (jo=0 and ATiin=0. The boundary conditions are written in the form of 
a matrix equation: MA = 0, where the components of the vector A are the amplitudes At2 and Ajiout, 
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respectively inside the PNR and reflected from it outside, and the 2x2 matrix M depends on two 
parameters - temperature, T, and PNR radius, /?. For example, for the j=l case, the QP amplitudes are 
obtained by solving the following coupled equations derived from Eqs.25a and b: 



(1) (1) 

A n -An 

Tlout TXout 72 72 

a (1) v'2 (-1) 



= at r = R 



d (1) 



^Tlout ^^ ^Tlout 



r ^Tlout"^ ^T2^^^T2 



n (-1) 



(26) 



Solving the determinant of Eq.(25) yields several solutions, but only the first ones, corresponding to the 
highest temperature, are plotted in Fig. 2. 




R{l.u.) 



Fig. 2. Temperature of the local instability for the case of Vortex, VTinst, and Quasi Polar, QPT,>,st , 
excitations, respectively, vs PNR radius /?. Total angular momentum j=l. Fl, F2 -fitting lines 
corresponding to y=125-1750/R^ and y =12234-324 A/R^, respectively. 

The temperature of instability is smaller for larger total angular momentum j> 1 and for a smaller PNR 
radius R. For example, for the case of a QP excitation, r,„st(/=2,/?=4)=73.8K and r,„st(/=2,R=6)=97.5K. This 
phenomenon is explained by the presence of a "centrifugal barrier", due to the angular momentum 
contribution ~j(j+l)/R^ in the Hamiltonian. The "average" momentum corresponding to a given density 
of displacement is larger for the case of a Vortex excitation than for the case of a QP one due to 
repelling from the PNR center. It is necessary to supercool the PNR below Td in order to compensate for 

the corresponding contribution to the free energy. This supercooling is stronger for the case of a Vortex 

r^ 



excitation, Ti^st ~ Tc2 - (1 ^ 0.25) 



a2*R^ 



. The respective displacements ^v and ^qp for V and QP 



excitations can be written in a simple form with spherical coordinates wheny=l, M=0, and at large r. 



^v 



exp(-kr) 
kr 



sin(6i)e^, f^p- 



exp(-kr) 
kr 



sin(0) eg, k = al^Tinst — Tc^/Cq 



(27) 



A Vortex (/=1) excitation describes the rotation of a rigid spherical shell around a polar axis and can be 
described as a manifold of closed lines (similar to the magnetic field lines in electrodynamics), which 
does not therefore possess an electric dipole moment. A quasi-polar excitation corresponds to 
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displacements parallel to the polar axis (17b) and integration over spherical angles reveals the existence 
of a non-zero electric dipole moment but only fory'=l. Following the previous Eqs (17b) and (26), the 
dipole moment, P, associated with the QP excitation can be written: 

(1) (1) (-1) 






Z 



M 



fSlo..™S2(''>-'2/3 + G^2 



(r)Vl/3]6'(i? - r) + 



;(1) 



;(-l) 



(28) 



GY{(R)[GY{^^^ir)^l2/3+G\.^'J^^ir)^'l/3]0ir-R)}, ^(x<0)=0, ^(x>0)=l 

where G's are defined in Eqs.AS, A7, and A8 in the Appendix. The unit vector x^ (H) defines the direction 
of the dipole moment P and is fixed extrinsically by interaction with chemical short range ordering, 
dislocations and so on. As shown in Fig. 3, the radial dependence of P^/ x,i is mostly confined inside the 
PNR and sharply decreases at the PNR surface. 

In Fig.3, we have calculated the density of TOW displacements for the V and QP excitations near 
Tins,. For QP excitations, we calculate separately the densities DQPpar for displacements parallel and 
DQPperp for displacements perpendicular to the radius r. These densities are defined as the square of the 
displacement averaged over angles inside a spherical layer. 



DV 



fj(kr) 



DQPpar 



G(^^(r) , DQPperp 



/-u 



(r) 



(29) 



the PNR surface looks as an almost impenetrable barrier for these displacements DQPpar parallel to the 
radius r is mostly concentrated inside the PNR such that as compare with the case of vortex excitation. 



T polarization 

DQP par 

— — DQP perp 
DV 




Fig.3 Radial layer distribution density of the normalized QP polarization, QP and Vortex displacements 
(30) near the temperature of the corresponding local phase transition, VTmspSO.l K, QPr,nst=112.23K. 
PNR radius /?=6,y=l. 

3.3 LOCALIZED AND QUASI LOCALIZED OPTIC MODES 

Local optic modes (LM) are characterized by an infinite lifetime as determined by lm(a))=0. They are in 
general discrete and exist in the "cage" created by the host media; the LM frequency is lower than the 
frequency of modes propagating outside the PNR, colm < Wo(l)- Therefore, LMs cannot be "directly" 
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excited by a propagating incident TO wave and only shallow LMs can scatter it effectively. Quasi 
localized modes, QLM, have an energy that is greater than the soft mode energy gap of the host 
medium. They result from the effective long-time interaction between an incident TOW and the PNR. 
They are seen as higher amplitude of the propagating TO mode in the vicinity of the PNR and can lead to 
its resonant scattering. The analysis of LM and QLM modes is presented successively in the following 
sections. 



Local optic modes 
Vortex LMs appear if a non-zero amplitude solution of Eqs (24a), (24b) exists when A ■,„ '°' =0. Similarly, 
QP LMs appear if a non-zero amplitude solution of Eqs (25a), (25b) exists when A Tiin=0. These become 
respectively: 

1 dgj _ 1 dgjou t 



dr 



1 da 






3j,out 



, r = R, lm{(ji)ij^) = 0, lm{qV) > 



dG. 



(1) 

Tlout 



4JU+^G. 



(-1) 
Tlout 



T2 



-(1) 
^Tlout 



dr 



^Tlout 



(30a) 
(30Z)) 



Two solutions of these, high and low frequency, are shown are shown in Figs 4a 



> 

0) 




T(K) 



Fig.4a. Frequency co of the Vortex and QP localized excitations, LM. The high (hf) and low (If) frequency 
branches are shown as a function of temperature. The PNR radius was chosen to be /?=12 andy'=l. The 
QP mode becomes unstable at T=121K. Curves Fl and F2 are fitting curves corresponding to y=(T-112)^^^ 
and y=(T-121)^^^ respectively. Also shown is the temperature dependence of the host medium soft mode 
gap, F3, which sets an upper limit for the value of the LM frequency. 

High and low frequency LMs correspond to different radial distributions of the displacement for a given 
value of the total angular momentum / (see later) and these frequencies increase with J. They are 
analogous to the different discrete levels corresponding to a given radial quantum number in quantum 
mechanics [17]. The temperature dependence of the low frequency LM branch, (OiMif, can be explained 
as follows. The typical momentum inside the PNR is proportional to b/R (b being a numerical constant) 
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and the LM excitation does not penetrate deeply into the host medium. Tal<ing into account expressions 



(9) and (5), we can write the frequency of the LIVl as : oiiMif ~ \{~Y + (<^2)^(T — TC2) or equivalently 

oHiMif ~ cl2^T — TC2, with Tc2 = TC2 — ( )^. The temperature dependence of the low frequency 

LM branch is thus similar to that of the soft-mode gap (^=0) but with a smaller value of Tc.. The 
frequency of the local QP mode is smaller than that of the Vortex for two reasons: 1) the QP wave field 
(/=1) contains a component with zero orbital momentum (15), and 2) QP mode interacts with the high 
energy LOW at the PNR surface. This interaction leads to a decrease of the QP mode energy, in 
accordance with the general results of second order perturbation theory. [17] Although not shown here, 
low and high frequency Vortex local modes correspond to different spatial distributions of the TOW 
displacement. Sample calculations for the following parameters (/?=12, j=l, 7=140 K, a)LMif=5.258, coLMhf = 
7.028, cjOo(1)=7.048, a)o(2) = 3.87 meV) reveal a maximum of the low frequency LM displacement density 
at r=6.6 and only one node at the PNR center. The high frequency Vortex LM displacement density has 2 
maxima at r=4 and r=11.2 and an additional density node at r=8.61 inside the PNR and the local h.f. 
Vortex mode is shallow in this case. This picture - increasing number of nodes with increasing energy - is 
similar to the dependence of wavefunctions on the radial quantum number in atomic physics [17]. 

The QP local mode structure is shown in Fig. 4b. 



1.0 -k^ 



0.6- 



■g 0-4. 



0.2 ■ 
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A DQPperp LM hf 
DQPpar LM hf 
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^0000000005, 



— I — 
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12 
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14 



16 



Fig.4b. Radial layer distribution density of the normalized QP displacements for the case of h.f. and l.f 
local modes. R=12,j=l, 7=140 K, QP CjOLMhf=5.123, (jOLMif=4.384 meV. 



Both DQPperp and DQPpar have a maximum at the PNR center. The displacement density for DQPperp 
of the h.f. mode has an additional node at r=6.41. The most unusual behavior is that of the h.f. and l.f. 
LM displacements parallel to r : both show no node and the density discontinuity at the PNR surface is 
stronger for the low frequency LMs. LMs exist within a limited temperature interval. For example, in the 
case of smaller PNR size /?=6, Vortex LMs only exist for T<140K. The structure of the TOW localized in 
PNRs could be studied by means of inelastic neutron scattering. 
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Quasi local optic modes 

It is impossible to directly "convert" a single propagating incident TOW to a local mode because the 
frequency of the incident wave is more than that of the LM frequency. However a TOW can be 
resonantly scattered by the PNR due to the existence of the so-called quasi local modes, QLM, the 
frequency of which appears in the continuous TOW energy spectrum. These modes are due to the long- 
time effective interaction between incident TOW and PNR, and lead to resonant TOW scattering. QLMs 
are similar to quasi-stationary states in quantum mechanics [17] and correspond to the solution of the 
dynamical equation, taking into account only the emitted (outgoing) waves Gcexp(iqir-i(jot). This solution 
is divergent at large distances from the center of the PNR, (lm(qi)<0), and its frequency o) contains a 
small negative imaginary part which corresponds to the QLM decay. QLM's are described by the 
solutions (24a,b), (25a,b) corresponding to Vortex and QP excitations respectively, jointly with the 
condition lm(cj)<0, lm(qi)<0. The QLM momentum, qi, and frequency are shown in Table 1 and Fig 5. 

Table 1. Parameters of a few Vortex and QP local (LM) and quasi local modes (QLM) for /?=6 and j=l. 



Excitation 


T,K 


Re(a)), 
meV 


-lm(a)), 
meV 


Re(ql), 
r.l.u. 


-lm(ql) 
), r.l.u. 


Re(q2), 
r.l.u. 


-Im(q2) 
, r.l.u. 


VLM 


130 


5.51 








-0.027 


0.08734 





VQLM 


150 


7.42 


0.668e- 
1 


0.190e- 
1 


0.715e- 
3 


0.783e- 
1 


0.174e- 
3 


VQLM 


170 


8.01 


0.387 


0.275e- 
1 


0.286e- 
2 


0.625e- 
1 


0.126e- 
2 


QPQLM 


130 


6.954 


0.026 


0.0239 


0.0107 


0.09407 


0.0027 


QPQLM 


150 


7.44 


0.0269 


0.024 


0.0126 


0.07875 


0.00385 


QPQLM 


170 


7.856 


0.0285 


0.02236 


0.016 


0.05864 


0.0061 
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Fig. 5. QP and Vortex local, LM, quasi local mode, QLM, frequency, and host medium soft mode gap u)o(l) 
as a function of temperature. R=6,j=l. 

Localized Vortex modes (V LM) exist for T<140 K and "convert" to QLMs for T>140K (Fig.5). Quasi polar 
LM and QLMs co-exist in the range 130K < T< 170K. 

Because the value of |lm(q2)/Re(q2)|is small. Vortex and QP QLMs correspond to "free" TOW motion 
inside the PNR. The parameters characterizing them strongly depend on temperature. 



3.4. TRANSVERSE OPTIC WAVE SCATTERING BY THE POLARIZED NANO REGION. 

The scattering cross-sections of transverse QP and V excitations have been calculated using Eq.(22). Figs. 
6- 9 illustrate their dependence on wavevector and temperature. 
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Fig. 6. QP and Vortex scattering by PNR (total and partial, j=l, cross section I) showing a strong 
resonance. The incident TOW momentum q^ is expressed in units n/a where o is a value similar to the 
lattice constant and /?=6, 7=130 K, Co=70 meV/r.l.u.; 1- marks the position of the QP quasi-localized 
state; the horizontal dash line marks the PNR geometrical cross-section. 
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a) 



b) 



Fig. 7 a). Total cross section of PNR as function of the incident momentum qi at the temperature T=130, 
150, and 170 K. a) for scattering by QP mode, long tail part of curves (qi >0.15) can be approximated by 
I(T=130K)=0.15/(7l^^ I(T=150K)=0.065 /qi^\ I(T=170K)=0.022/qi^^ b) for scattering by for scattering 
by Vortex mode, the long tail part of curves (qi>0.15) can be approximated by I(T=130K)=0.15/qi^^, 
I(T=150K)=0.075 /ql^^ I(T=170K)=0.02/(7i^^ 
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Fig. 8. Vortex partial-wave scattering cross sections, corresponding toy=l, 2, 3, 4 and R=6 a) atT=130 K 
and b) 95K, We note in particular the dominant role of the quadruple (1=2) scattering. 
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Fig. 9. Partial (/=1) QP and Vortex scattering cross section as function of qi with angular momentum at 
the temperature 130, 150, and 170 K. R=6, Co=70 meV/r.l.u.; 1,2, and 3 - mark the positions of the QP 
and Vortex quasi-localized states. 

The TOW scattering obviously also depends on the PNR radius (Fig.lO). 
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Fig.lO. QP and Vortex total scattering cross section as a function of PNR radius value R : ql=0.1; 7=130 
and 190K. 

Diffraction maxima and minima are seen here for the both V and QP scattering at T=130K. Such maxima 
and minima are absent at T=190K because the soft-mode gaps inside and outside the PNR have almost 
of the same value at higher temperature (see Fig.l) and the scattering there is weak. The moderate 
oscillations in the curve corresponding to Vortex scattering at T=130K are created by partial harmonic 
contribution to the cross-section. 



3.5. DISCUSSION 



The calculated results shown in Figs. 5 - 9 can be interpreted qualitatively as follows. 

For the case of very small incident momentum qi, the TOW dynamics is similar to the scattering of a 

slow particle in quantum mechanics [17]. The calculated value of the scattering cross-section is very 

small, a-q/for q':[<0.01 andy'=l, which is in agreement with the general theory of scattering [17]. 

The scattering cross section I decreases at large momentum, qi>0.15 (Fig. 7a, b). For the case of QP and 



Vortex excitation I(T=130K) = 0.15/qi ^ , I(T=150K) = 0.055 /q^ , I(T=170K)= ) = O.OII/qi and 
I(T=130K) ) = 0.15/(7i^ ^ I(T=150K) ) = 0.055 /qi^ ^ I(T=170K) ) = 0.022/qi^ ^ respectively. 
Therefore the value of the potential scattering cross section has a wide maximum in the intermediate 
region around qi ~ 0.05 ■;- 0.1. The reason for the decrease of the cross section at large momentum is as 
follows. Scattering is due to the difference in the soft mode gaps (Eq.5) inside and outside the PNR, 21a» = 
u)o(l) - u)o(2), which is equal to 4.57, 2.3 and 1.08 meV at 130K, 150K and 170 K respectively, and defined 
by the value of the dimensionless parameter /icj/Coqi <1. This potential scattering is stronger at low 
temperature and decreasing for large values of q^. The effect of scattering disappears at the 
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temperature Ttransp- (see Figl) at which the soft-mode gaps inside and outside the PNR become the same 
and the soft-mode is no longer scattered by the PNR, Ttransp~194 K Then, QP and V scattering is small 
and almost the same at 7=190K (Fig.lO). 

However this picture of the potential scattering does not exhaust the problem. As shown in Figs. 6-9, 
strong resonant scattering of the TOW is clearly seen at small q. There, the value of the scattering cross 
section is 4 ■;■ 8 times larger than the value of the geometrical PNR cross section nR^ -113 \.\i}. 
The difference between QP and Vortex scattering is essential in the resonance region, due to the 
effective long time interaction between excitation and PNR. QP and V cross sections are similar outside 
of this region where the corresponding cross section can be calculated in the Born approximation and 
with a plane wave basis. 

The calculated total scattering cross section as function of the momentum qj exhibits a strong peak 
asymmetry and bumps. These peculiarities are the results of scattering QP and V excitations with 
different angular momentum j (see especially Fig. 8). These bumps are more visible at low temperature. 
The relative value of the partial cross section depends on temperature. For example. Vortex peaks at are 
smaller than QP ones T=130 K and 170 K but not at T=150K case (compare Figs. 7a and b). This is due to 
the presence of quasi local (QLM) and local (LM) modes. As explained above, the propagating incident 
TOW can "directly" excite QLMs, leading to strong resonance scattering. Similar resonant scattering is 
not possible by LMs because their energy is smaller than the energy of the incident TOW. In general, 
only shallow LM will have a significant effect on scattering. The difference in the effect of LMs and QLMs 
on scattering is confirmed by a comparison of scattering corresponding toy=l (Fig.9) and data shown in 
Fig. 5 calculated for the same casey=l. QP LMs are not shallow and probably do not result in a strong 
scattering effect. The ql position of the QP quasi local modes position is almost independent of 
temperature and is found near the scattering peak at T=130 and 150 K (Fig. 9). By contrast, the position 
of the Vortex QLM is found at small ql, Re(ql)=0.019, and near the strong scattering peak maximum at 
T=150 K, but it moves away at higher temperature (Re(ql)=0.0275 at T=170 K) when the role of 
resonance scattering and the value of the cross section correspondingly decrease. Only shallow Vortex 
LMs exist at low temperature, T=130K, leading to only moderate V scattering. 

Finally, we should note that the position of the V cross section maximum along the ql axis depends on 
temperature unlike that of the QP cross section (Fig.9). This fact is confirmed by the position of the QP 
and V local and quasi local modes in Fig.5. 

In Figs. 5 and 9, the fact that the QLM position and maximum of the scattering cross section do not 
coincide is due to interference between resonance and potential scattering. In Fig.8a, the partial cross 
section peaks corresponding to angular momentum />1 appear at higher qi positions because the 
centrifugal force repels the TOW from the PNR. In practice, resonance scattering is most important for 
7=1 and y=2. The appearance of QLM for j=2 leads to the dominant role of the Vortex quadrupole 
scattering at 7=96K (Fig.8b). 

4. CONCLUSION 

We have developed a model describing the dynamics of the Transverse Optic Wave (TOW) in a medium 
containing one spherically symmetric Polarized Nano Region (PNR). Our analysis is limited here to the 
temperature region below Burns temperature but above the temperature when PNRs become static or 
quasi-static polarized nano domains (PND). Two separate types of excitations are found. Vortex and 
Quasi-Polar (QP). Vortex excitations correspond to closed polarization lines while QP excitations 
correspond to open polarization lines with a dipole moment (if the total angular momentum j=l). 
Standard boundary conditions are applied at the PNR surface in the case of Vortex excitations. For quasi 
polar excitations, after excluding the Longitudinal Optic Wave (LOW) contribution, the boundary 
conditions correspond to continuous TOW displacements but discontinuous derivatives. The 
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discontinuity is connected with near field effects at the PNR surface where the transverse optic mode 
mixes with the longitudinal one. 

We show that a local violation of stability ("phase transition in PNR") occurs at a temperature Tipst, the 
value of which is smaller in the case of Vortex than QP excitations. Therefore, the QP condensation 
occurs at a higher temperature than the Vortex condensation upon cooling. We also find that a localized 
TO mode (LM) can exist inside PNR and calculate its LM frequency, cjlm- A local phase transition then 
takes place when WtM=0. We have computed the cross-section for scattering of a TOW by a PNR and 
found a strong resonance. Vortex and QP scatterings are practically the same outside the resonance 
region. In the resonance region however, the Vortex scattering is significantly different from the QP one, 
which is attributed to the presence of localized and especially quasi-localized modes. The structure of 
these modes could be studied by means of inelastic neutron scattering. 

A detailed discussion of the experimental situation is out of scope for our paper. We note that 
Gehring et al [6, 9] found evidence of a strong step-like increase in TOW damping (waterfall) for an 
incident TOW momentum qi<qo- The value of 1/qo was found to be consistent with the size of PNRs. We 
also find a broad TOW scattering peak around momentum qo- In addition however, our calculations 
reveal the presence of a strong increase in TOW scattering at small q^, due to the presence of shallow 
localized and especially quasi localized states within and around PNRs. The value of the calculated 
resonance scattering cross section is 4 ■;- 8 times greater than the PNR geometrical cross section. At small 
qi the model of scattering by individual and well separated PNRs is no longer valid since interference 
between waves scattered by different PNRs becomes essential. Our calculation shows that strong soft 
mode scattering is present in the paraelectric phase. Probably, in order to have a waterfall, it is not 
necessary to assume that PNRs are in a condensed, ferroelectric phase (PND); it is sufficient that the 
system contain micro regions, each of size " 10 -^ 20 lattice constants, with values of the soft-mode 
energy gap different from that of the host medium. This, of course, does not exclude the important role 
of PNDs in TOW scattering at lower temperatures. However, our calculations are limited here to the 
temperature range 7'>T|nsf An analysis of the "condensed phase" (7'<Tinst) is more complicated since 
anharmonic interactions will be important in that case. 

APPENDIX. 

DEDUCTION OF BOUNDARY CONDITIONS FOR THE CASE OF QUASI POLAR EXCITATIONS 

First of all we discuss boundary condition terms (7b), (7c) including the electrostatic potential <p. We 
have solution of Eqn. (7b) outside of PNR surface taking into account Eq-s (19, 20a, 20b) 

* 

8 2 dSGjir) 
^ ir ^ ) - j{j + \)SG. (r) = (Ala) 

Solutions of (Ala) which is limited inside and outside PNR are well known 
SGj{r)^J'PrJe{R-r)+A''P_^r~J~^e{r-R),e{x<Q)^Q, ^(x>0)=l (Alb) 

where A s^\ A s,-i^' -arbitrary constants and R- PNR radius. 
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The boundary condition parts of (7b), (7c) lead to Eqs (Ale), (Aid) 

(2) j-\ * (2) (1) -j-2 

A jR - Ane f. , ^ -A (j + \)R 

o J,-l o,-\ 

* (1) * (1) 

- Am f . . , - Ane f . , (Ale) 

] ,in ,-\ J, out, -I 

* (2) (2) j * (1) (1) -j-\ 

4716 G . (R) + A^ R ^ Am G . (R) + A ^ R (Aid) 

JO J d,-\ 

where /^'j^.i, f^'jjn^.i, f^'j^out,-! - transverse optic wave displacement component along the PNR radius. It is 
interesting to note that longitudinal optic wave amplitudes are absent in (Ale). We can solve (Ale, Aid) 
with respect to A J^' and A s,-i^' and satisfiy the BC's (7c) and (7d). The electrostatic potential (p is 
included in Eqn. (7a) in the form d(p/dra and therefore does no affect on the circulation curl(^ dynamics 
at the TOW. Due to the equality A6(p=0, the electrostatic potential component, 6(p, does no affect on 
the electric charge e*cliv(^) dynamics at the LOW. Therefore we omit boundary conditions (7b), (7c) for 
the electrostatic potential in the reasonably approximation. 

In order to satisfy other boundary conditions, it is necessary to take into consideration functions (17a, b, 
c), (19). These functions contain two components, proportional to the function y'^'jm and Y^'^'jm 
respectively. The contribution to the optic wave displacement of the corresponding components can be 
written as follows: 

^ = G (r)Yj^ + G (r)yjM (^2) 

The "force" terms of the boundary conditions (7a) contain expressions (A3) and (A4). 



(-1) ^/7(7^T) (1) d (-1) 2 (-1) 
'a " JM ,a 



^'■^(^)«. = ^,-A^.[ r G +^G +-G ] (A3) 



(1) SG^^) v7(7^ (-1) 1 (1) 

(-1) 5G^-1) V7(7+l) (1) 2 (-1) 

The boundary conditions (7a) can be formulated as a continuity requirement on the value Oaene 



(1) 2 8G''\^UU+ll (-1) 1(1). 
Y^-;P {c2[^^-^/^^Mc,(l) + 2^(-l)]+c2[2VZMc,(l)_4c,(-l)]| (A5) 

jM ,CiC L Qf f f If f 

For the transverse quasi polar (17b) and longitudinal (21) optic modes, the corresponding expressions 
are respectively: 
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= G r = ^ .^^ ["^7+1 ^^'"^ + ^7-1 ^^''^^ 



G 



G(l).G(l)^.G^.^L7i7±^, G^-i;.G^-U 



L = 



dr 



(A6) 

(A7) 



The functions gj(l<r) are described in (16). The value of an ion plasma frequency (8b) is very large, Qp 



:(-i) 



»cjo , Qp »(jOo(1,2) . Therefore \G^ \ « \G^ ''\ and the LOW displacement is directed along the PNR 
radius when Qp->°°. The functions Gj, Gyp, Gjout can be written as follows: 

3/2 j 
r < R, G^. = (2;r) / ^y+1/2 (^^ (2)r) / V^^ (2)r, 

r>R, Gjjn ={2k)^'^ '^^j+1/2 ^^^ {2)r)l ^ Ki{2)r , 
r>R, Gy,oM?=(2;r)3/2.7^(U^^(^^(2)r)/V^^(2)r, 



K (\,2) = /^/f3^+ft;Qn.2;-»^/C 



Z 



« /Q / C 



Z, 



Kl 



R » 1 



(A8) 



Finally, the full boundary conditions (A9a, b, c, d) contain the optic wave amplitudes At2, Al2 (inside PNR) 
and the amplitudes Ajiin, Ajiout, Aliqui (outside PNR). The momentum Ki(l) is imaginary and we therefore 
have ALiin=0. 
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We are interested by the case when a value of LOW frequency Qp^°° and |Kl.RI»1. Let us suppose 
that Al2Gl2'^', ALioutGiiout'^' "I/Kl so that we can omit terms "" Al2, Aaout in (A8a) (but not in (A8b, A8c)). 
We find expression (AlO) from (A9b), (A9c) and (A7). 
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We can omit terms ~ Aqoui , l^n, in (AlO) and find the boundary conditions (26a), (26b) used in the main 
text of the paper. 

The contribution of the LOW to the displacement parallel to PNR surface is small, "l/ 1 Kl| . However the 
contribution of the LOW to the displacement perpendicular to PNR surface is essential within a narrow 
layer of thickness "I/IKlI. This component of the LOW displacement leads to the existence of a non- 
zero term on the right hand side of Eq (26b). 

In order to check the quality of our approximation |Kl|R^ °° we calculated the QP TOW scattering 
cross-section a exactly (Eqs. A9a, b, c, and d) for the case Qp=30 meV, /?=6, 7=130 K, Co=70meV/r.l.u., 
Q=105 meV/r.l.u. We found that this exact value of a is in good agreement with the value of I (Fig. 5) 
calculated in the approximation | Kl| R^ °°. 
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